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Abstract —Chirped sinosoids and interferometric phase plots 
are functions that are not periodic, but are the composition 
of a smooth function and a periodic function. These functions 
functions factor into a pair of maps: from their domain to a 
circle, and from a circle to their codomain. One can easily imagine 
replacing the circle with other phase spaces to obtain a general 
quasiperiodic function. This paper shows that under appropriate 
restrictions, each quasiperiodic function has a unique universal 
factorization. Quasiperiodic functions can therefore be classified 
based on their phase space and the phase function mapping into 
it. 

I. Introduction 

Chirped sinosoids and interferometric phase plots are func¬ 
tions that are not periodic, but are the composition of a 
smooth function and a periodic function. These functions 
functions factor into a pair of maps: from their domain to 
a circle, and from a circle to their codomain. Functions that 
are the composition of a smooth map and a periodic function 
can be rather different from almost-periodic functions, since 
chirps do not lie in the closure of the space of trigonometric 
polynomials with any reasonable norm. Additionally, chirps 
are not periodic to first order, so they are not amenable to 
asymptotic approaches like the famed WKBJ method H). 

One can easily imagine replacing the circle with other phase 
spaces. This paper shows that under appropriate restrictions, 
each such quasiperiodic function has a universal factorization. 
These quasiperiodic functions can then be classified canoni¬ 
cally based on their phase space. We expect that this canonical 
representation could serve as the starting point for practical 
filtering algorithms for quasiperiodic signals. 

A. Contributions 

This paper presents a general definition of a quasiperiodic 
function and factorizations of it that ensure the phase space 
is a manifold. Given this definition, it proves that unique 
universal factorizations of quasiperiodic functions exist. To 
connect the theory of quasiperiodicity to periodic functions, 
it completely characterizes universal S' 1 -quasiperiodic factor¬ 
izations of periodic functions. This paper ends by presenting 
an initial algorithmic framework for producing factorizations 
of quasiperiodic functions. 

B. Historical context 

The problem of recovering the phase of an analytic signal 
from its amplitude is called “phase retrieval”, and has occupied 


the attention of researchers at least since the 1950s in the 
context of x-ray crystallography & Work continues even 
recently (for instance |[3|), though most of the results relevant 
to this article were obtained in the 1980s 0- Numerous 
methods, such as those based on iteration [5J, asymptotics 0, 
or gradient search methods 0, have been developed to address 
the problem, though they all rely heavily on the assumption 
that the phase lies on a circle. Although phase ambiguities are 
prevalent in time series data, they are quite rare in the case of 
a two-dimensional domain 0. 

The study of non-circular phase spaces has been pursued 
most aggressively by researchers in dynamical systems. Most 
of their approaches are based on the famous embedding 
theorems of Whitney 0. In particular, the Takens delay¬ 
embedding theorem (see for instance El, ED) can obtain 
topologically-accurate models of phase spaces that are suitable 
for factorizations of quasiperiodic functions. 

Our approach is inspired by the work by DeSilva et al. 
ED, which uses persistent cohomology to construct circular 
coordinates. They use a delay embedding approach to obtain 
a model of the underlying phase space, followed by persistent 
cohomology to identify nontrivial 1-cycles, which they smooth 
using harmonic 1-forms. The class of functions they consider 
are precisely those that factor into a diffeomorphism and a 
periodic function. We generalize beyond their framework in a 
complementary way, by exploring more general factorizations 
and justify that their delay-embedding approach still works in 
our case (though not necessarily the harmonic smoothing). 

II. Motivating examples 

Every smooth periodic function / : R —> TV from the real 
line to a manifold N can be written as the composition Fo(f>, 
where tj) : R — > S 1 is a covering map. The function F : 
S 1 -A- N encodes the behavior of / on a fundamental domain. 
The function tj) extends F periodically from .S' 1 to the whole 
real line. Since tf> is a covering map, it is a surjection and a 
local homeomorphism. Thinking about the topological spaces 
involved, S 1 can be written as the quotient R/<j>, which we 
call the phase space for the factorization F o tf>. 

There is, however, a different quotient that may be taken as 
well, namely R//. This quotient serves as the phase space for 
a different factorization f = G on, where n : R —>■ R// is 
the canonical projection. In this factorization, n is surjective 
and G is injective. 


Since f = F o f, any pair of points x, y £ R with n(x) = 
n(y) will also satisfy f(x) = f(y). Therefore, R// might 
appear to be a better object to work with than R/f - since 
it is the minimal specification of the function / - but R// 
isn’t usually a manifold. Wherever / has a critical point, the 
corresponding point in R// will usually not have a Euclidean 
neighborhood. 

Example 1. Consider the function / : M — >- R given by 
f(x) = sina:. Then R// is the closed interval [—7r/2,7 t/ 2], 
which is not a manifold. (It is a manifold with boundary, 
however.) 

III. Quasiperiodicity 

In order to ensure that the phase space of a factorization 
of a periodic function is a manifold, we need to constrain 
the projection to the phase space. The lesson from differential 
topology is that this projection should be a surjective submer¬ 
sion (see |12j Prop 7.16-7.19] for instance). 

Definition 2. Suppose M, N, C are smooth manifolds. A 
smooth function u : M —> N is called C-quasiperiodic if it 
has a factorization u = Uo<j> where <f> : M —> C is a surjective 
submersion. Such a factorization is called a C-quasiperiodic 
factorization. We call C the phase space of the factorization. 

Quasiperiodic factorizations of a given function are not 
unique, even for a given phase space. 

Example 3. Consider the periodic function u(x) = sin 2nx as 
a function M —>• R. One quasiperiodic factorization uses R as 
the phase space and </> as the identity function. The function 
u also has many .S' 1 -quasiperiodic factorizations. Using the 
phase space C = S l = R/Z, one factorization has <j>(x) = 
cc/2 mod 1 and U{9) = sin(47r0). Another factorization is 
4>'{x) = x/3 mod 1 followed by U(9) = sin(67rd). 

Fortunately, there is a canonical notion of a quasiperiodic 
factorization, namely one that satisfies a universal property. 
This neatly balances the need for a minimal specification with 
the requirement that the phase space be a manifold. 

Definition 4. A C-quasiperiodic factorization u = U o f of 
a smooth function M —> N is universal if for any other 
quasiperiodic factorization u = U' o ft with </>' : M —> C ', 
there exists a smooth function c : C' —> C such that the 
diagram below commutes: 


M-^C 



C-u* N 


A universal quasiperiodic factorization has the minimal 
phase space among quasiperiodic factorizations. It therefore 
generalizes the notion of a fundamental domain for a periodic 
function. This connection is made precise in Theorem [8] 


Theorem 5. Every quasiperiodic function has a universal 
quasiperiodic factorization, which is unique up to a diffeo- 
morphism applied to the phase space. 

The proof of this theorem is given in Section [VJ and a 
complete characterization of universal S' 1 -quasiperiodic fac¬ 
torizations of periodic functions is given in Section IIVI 

Remark 6. For a function u : M — > N, the quotient M/u 
also satisfies a universal property like the one described in 
Definition [4] However, as noted in the previous section, M/u 
is generally not a manifold. 

Because of their uniqueness, universal quasiperiodic factor¬ 
izations govern the structure of all quasiperiodic factorizations. 
Therefore if we find a quasiperiodic factorization by any 
means, we know its structure is constrained from the outset: 
it must have the universal phase space as a submanifold and 
must descend to this submanifold as a quotient. 

In some cases, universality of a quasiperiodic factorization 
is easy to detect as the next Proposition shows. 

Proposition 7. Suppose u : M —> N is a constant function. 
Then u : M —> * —> N is its universal quasiperiodic 
factorization, where * is the space with one point. On the other 
hand, if v = V o cf> is a quasiperiodic factorization where V 
is injective, then it is a universal quasiperiodic factorization. 

Proof: For the first statement, consider any other factor¬ 
ization u : M —> C —> N. Construct the unique function 
c : C —> * and observe that the diagram shown below left 
commutes: 


M -^ * M —^ C 



C -^ N C ->■ N 

V' 


For the second statement, consider the diagram above right 
and suppose that v = V' o ft is another quasiperiodic 
factorization. Since V is injective, it has a partially defined 
inverse U _1 : v(M) —» C, which means that c : C —>• C 
given by c = V ~ 1 o V' completes the commutative diagram. 

■ 

IV. Universality for periodic functions 

This section presents a complete characterization of the 
universal S' 1 -quasiperiodic factorizations of periodic functions, 
thereby establishing a link between the theory of quasiperiodic 
functions and the traditional theory of periodic functions. The 
phase space of such a universal quasiperiodic factorization is 
precisely its fundamental domain. 

Theorem 8. Suppose u : R —> N is a smooth, non-constant 
T-periodic function. Then 

1) u is S 1 -quasiperiodic and has a phase function (j> : R —> 
S 1 = R/Z given by f(x) = x/T mod 1. 

2) The factorization is universal unless there is an S 1 - 
quasiperiodic factorization of u whose phase function 
is non-injective on an interval of the form [a, a + T). 








Example 9. As an example of a universal quasiperiodic fac¬ 
torization of a periodic function, consider the S' 1 -quasiperiodic 
function defined on [0, 2] given by the factorization f = Uo<j> 
where 

f(x) = x/2 mod 1 

and 

sin47T0 ifO<0<l/2 
i sin SttO if 1/2 < 0 < 1. 

This is an example of a universal quasiperiodic factorization. 
Notice that any factorization that has a phase function with a 
period smaller than 2 would be unable to correctly reproduce 
the change in amplitude of u. 

Example 10. The smoothness of u is a necessary condition 
for Theorem [H] Consider the case of the function given by the 
factorization g = U o cf> where 

cj)(x) = x/2 mod 1 

and 

sin47rd if 0 < 9 < 1/2 
sin87rd if 1/2 < 6 < 1. 


1) Observe that [0. T) is a fundamental domain for u. 
Also note that cf>(x) = x/T mod 1 is a surjective 
submersion. Hence with U{9) = u(9/T) for 0 < 9 < 1, 
then u=U o cf>. 

2) Suppose that u has another quasiperiodic factorization 
u = U' o ft so that ft : R —»• C is a surjective 
submersion. Because u is not constant and ft is a 
submersion, rank ft = 1. Therefore ft has no critical 
points. Since ft is surjective, we observe that C must be 
a manifold of dimension 1 without boundary. Therefore, 
C is of two possible types: 

a) C is diffeomorphic to (a, b) C R. Because qY 
has no critical points, it must be monotonic and 
therefore injective. Since it is a submersion, qY 
has a continuous inverse by the inverse function 
theorem. Construct a map c : (a, b) —> S 1 by 
c(x) = ft~ 1 (x)/T mod 1. This map is the re¬ 
quired one, since 

c{ft(y)) = (<?%))/ T modi 

= y/T mod 1 = <f(y). 


almost as in the previous example, though g is not smooth. 
However, there is also a new factorization g = U' oft 

x if 0 < x mod 2 < 1 

2x — 2 if 1 < x mod 2 < 2, 

with 

U' (9) = sin27rd. 


b) C is diffeomorphic to S 1 . In this case, we can 
simply write <// : R — > S 1 . By hypothesis, </>' is 
injective on every interval of the form [a, a + T). 
Since u = U' o </>' is periodic with period T, </>' 
must also be periodic with period a multiple of T. 
Thus we simply define c : S 1 —>■ S 1 by c(x) = 
4>'~ 1 (x)/T mod 1. 


Notice that </>' has fundamental period 2, but that there is 
no map c that makes the appropriate diagram commute - 
contradicting the universality of g = U o <fi. 

Proposition 11. Any function that is the composition of a 
diffeomorphism D : R —>■ R and a periodic function p : R —>- 
N is S 1 -quasiperiodic. 

Proof: The proof is a simple commutative diagram: 


p 



(fi'=(f)oD 


Lest one think that S' 1 -quasiperiodic functions are all es¬ 
sentially just periodic functions with distorted timescales, the 
converse is not true. 

Example 12. Consider the function it : R->- S' 1 = R/Z 
given by u(x) = arctanx mod 1. This function is S 1 - 
quasiperiodic - trivially - since arctan x is a submersion and 
has a range larger than [0,1). This is completely different than 
a periodic function since a preimage of u contains at most 3 
points - not an infinite set. 

Proof: (of Theorem [8} 


V. Existence of universal quasiperiodic 

FACTORIZATIONS 

Lemma 13. If u = U o <fi = U' of are two universal 
quasiperiodic factorizations with phase spaces C and C, 
respectively, then C and C are diffeomorphic. 

Proof: Observe that because the factorizations are both 
universal, there are smooth maps c : C —> C' and d : C' —> C 
making the diagrams below commute. 


( f> (h 

M C M C 



C -^ N C' -^ N 

u’ u' 


because both <p and ft are surjective, this implies that both c 
and c' are surjective. 

Suppose x,y £ C such that c(x) = c(y) G C'. Since f 
is surjective, there are x',y' G M so that </>(x') = x and 
4>{y') = y. Since the diagrams above commute, d oft = ft 
which means that c'(ft(x')) = x and c'{ft(y')) = y. However, 
ft{x') = c(x) = c(y) = ft{y'), so we must conclude that 
x = y. Thus both c and c! are bijective and are each others’ in¬ 
verses. Since they are both smooth, they are diffeomorphisms. 














Lemma 14. Suppose that u = U o <fi = U' o 0' are both 
quasiperiodic factorizations of the same function u : M —>■ 
N with phase spaces C, C' respectively. Then there is a 
quasiperiodic factorization u = U" o f" and maps p, p' such 
that the following diagram commutes: 


C 



Proof: Consider C" = (CUC , )/~ where x £ C ~ y £ 
C' if there is a z £ M such that <j>(z) = x and f(z) = y. Let 
p : C —>• C" and p' : C' —> C" be the canonical projections. 
Then the function cf>" : M —»• C" is given by f" = p o = 
p' o (f. It is well-defined by construction: if z £ M, then 
</>(z) ~ 4>'{z) which have the same images in C" through p 
and //. 

<f>" is surjective because both <fi and ft are, since if [x] £ C" 
for some x £ C, then surjectivity of </> implies that there is a 
z £ M with 4>(z) = x. 

What remains to be shown is that C" is a manifold and 
that ft' is a submersion. To this end, let [x] £ C". Without 
loss of generality, assume that x, y £ [x] for some x £ C and 
y £ C'. Let z £ c/r 1 {x). By the constant rank theorem ( IT2l 
Thm 7.13]), the tangent plane at z splits 

T Z M = {d z ft)~ l {T x C) ®K = {d.fty 1 {T y C') © K\ 

where K, K' lie in the kernel of the derivatives of </> and ft. 
Therefore, only vectors in [d z ft)~ l (T X C) Cl ( d z ft)~ 1 (T y C') 
correspond to the tangent plane to C" at [x] - hence C" 
is a manifold. Any other vector in T Z M corresponds to a 
direction that does not displace either x £ C or its corre¬ 
sponding y £ C'. Such points are all attached to [x] £ C". 
This also explains that p is a surjective map from T X C to 
[d z ft)~ l {T X C) Cl ( d z ft)~ l (T y C'). Hence ft' is a submersion. 

■ 

Proof: (of Theorem [5]) Suppose that u : M —> N is a 
quasiperiodic function. Let I be the set of all quasiperiodic 
factorizations of u so that for each i £ /, u = Ui o fa with 
<j>i : M —» Ci and U t : Cj —> N. This set has a partial order 
< in which i < j if there is a map c l3 : Ci —> Cj such that 
the diagram below commutes 



Notice that < is transitive because if i < j and j < k, then 
Cifc = Cjk o dj, so i < k. In this way, we construct the phase 
space for the universal factorization by 

C* = lim Ci = |_| Ci/ ~ 

ie/ 


where x, £ Ci ~ Xj £ Cj if there is some k £ I such that 

C-ikipti) — ('ttj ) • 

The function (f> t : M —> C* can be constructed by taking 
the composition of any 0 , with the projection Cj ^ C. t , and 
the function 17* : C* —» is given by the constraint that 
u = 17* o (j} t . Notice that 17* is well-defined by construction, 
since points in Ci and Cj are identified only when they have 
the same image in N. 

Finally, we note that the proof of Lemmafl4leasilv lifts from 
a pair of quasiperiodic factorizations to an arbitrary collection 
of factorizations. Thus Cj is a manifold and </>* is a surjective 
submersion. ■ 

VI. Practical factorization algorithms 

For S' 1 -quasiperiodic functions that arise from the com¬ 
position of a diffeomorphism and a periodic function on R 
(not the case of Example IT2T) . DeSilva and his collaborators 
present an effective factorization algorithm using persistent co¬ 
homology CD. Their algorithm is based on a delay-coordinate 
embedding approach (see 0 , m for instance). Suppose that 
u : R ra —> N is a quasiperiodic function for which we 
would like to obtain a quasiperiodic factorization. Select a 
collection of k vectors, {iq,..., Vk} and define the function 
f : R 71 —>■ N k+1 by 

</>(x) = (u(x), u{x + vi),..., u(x + v k ))- 

and the function U : N k+1 —> N by U(xo,..., Xfc) = Xo- 
Clearly, this is a factorization u = U o 0. According to a 
generalization of Takens’ theorem iflOl Thm 2.7 and Rem 2.9], 
a generic collection of {vi} yields a cj> which is an immersion 
if k > 2 n. If we restrict consideration to <t>(R n ) C N k+1 , then 
usually f is a surjective submersion. Therefore, this approach 
yields a quasiperiodic factorization of u. 

Although this approach yields quasiperiodic factorizations, 
they are often quite far from being universal. For instance, 
consider the the S' 1 -quasiperiodic function u(t) = sin(l/f) 
for t > 0. With k > 3, the delay-embedding approach will 
yield a quasiperiodic factorization, but this factorization will 
have R for a phase space, not S 1 as it should be if it were 
universal. 

Other approaches for identifying phase spaces have been 
suggested as well, such as the use of successive derivatives 
(u(x), ^j(x), ^t(x), ... ^ in lfT3l . but these too tend to 
recover R as a phase space. 

VII. Conclusion 

This paper proposes a new definition for quasiperiodic 
factorizations of functions and proves universal factorizations 
exist. It characterizes these universal factorizations in the case 
of periodic functions. Finally, it offers an algorithmic approach 
for producing quasiperiodic factorizations, which may fail to 
be universal. It remains to discover factorization algorithms 
that ensure universality. 
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